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, For the first time two separate production routes of Tc radionuclides have been studied bom- 

barding 7 Li on nat Zr and 9 Be on 89 Y. Excitation functions of the evaporation residues produced 
' in those reactions have been measured using stacked-foil technique followed by the 7-spectrometric 

D ■ studies in the energy range 37-45 MeV and 30-48 MeV respectively. Measured excitation functions 

have been compared with those calculated using the nuclear reaction model codes PACE-II and 
ALICE91. Experimental results show good agreement with the theoretical predictions. Compound 
nuclear reaction is the key mechanism in producing evaporation residues. 
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^ ' I. INTRODUCTION 

a/ 

Tc is the element of interest since last three decades because of its several elegant features 
and practical applications e.g., in the field of nuclear. Out of identified isotopes of Tc, 21 are 
On : proton rich, having (3 + and/or e decay mode, while rest 20 are neutron rich, mostly (3 emitter. 
Neutron rich Tc isotopes are short-lived except 98 Tc (T 1/2 =4.2 My) and "Tc (T 1/2 =2.111 x 10 5 
y). The only odd-even radioisomer, 99m Tc is extensively used in nuclear medicine because of 
its half-life (T 1 / 2 =6.01 h) and 140.474 keV 7-ray, both suitable for in vivo imaging. In fact, 
every year millions of people throughout the world are diagnosed using 99m Tc. However, neutron 
rich Tc radionuclides, even 99m Tc are not well suited because of extremities in terms of half 
lives for laboratory experiments of tens of hours duration. On the other hand proton rich 7 
and P+ emitting Tc tracers ( 93 Tc, 9 4™, 9 4 TC; 95m,95 TCj 96 Tc ) exn ibit required salient features to 
serve diverse applications, particularly in nuclear medicine. For example, 94m Tc (Ti/2 = 52.0 m, 
I / 3+=70.2%, = 2.44 MeV) is potentially promising radionuclide for PET (Positron Emission 
Tomography) imaging of organs because of its high positron emission energy. Short lived 7-emitters 
93 Tc (T 1 / 2 =2.75 h) and 94 Tc (T 1 / 2 =293 m) are convenient for short span experiments while 7- 
emitter 95 Tc, which has comparatively long half life of 20 h, can be useful to keep track of the 
metabolic function of human brain and heart. Therefore, information on the nuclear data for the 
production of proton rich Tc tracers having considerable half life is important. 

Several studies have already been reported which aimed to enquire the nuclear properties of 
proton rich Tc isotopes as well as isomers. Usually proton rich tracers are produced in medium 
energy accelerators usinglight ion (p, d, 3 He, a) induced reactions on molybdenum (Mo) or niobium 
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is mononuclidic in nature. Therefore, investigation of nuclear properties of Tc isotopes, which are 

Soduced by a and 3 He induced reactions on 93 Nb, has got importance in basic nuclear physics 
Large number of reports are available on the production cross section of 
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Tc radionuclides from proton induced reactions on enriched isotopes of Mo covering wide energy 
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range 
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IS III lilt- Excitation functions f 93 - 93m ' 94 > 94m ' 95 ' 95m > 96 Tc from 



proton and deuteron induced reactions on natural Mo have been reported in ref [19ll20l.l21l.l22l.l23|. 

Moreover, cross sections of 95 ' 96m ' 96 Tc have also been measured by neutron induced reactions on 
enriched 96 Ru as well as natural ruthenium target 
and 7-ray induced reactions on 99 Tc have been reported in 



271 ]. Even production of 95 > 96 Tc by proton 



2b. 
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However, cross section data of proton rich Tc isotopes produced via heavy ion induced reactions 
is rare. A few literatures are available where excitation functions of Tc isotopes are measured from 
12 C induced reactions on 93 Nb and 
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3l| aiming to study the reaction mechanisms involved 



in the particular reaction. Unlike neutron and light charged particle reactions, understanding of 
reaction mechanisms as well as nuclear data involved in the heavy ion reactions are not adequate 
till date and hence demands much attention in this direction. 

In order to satisfy the growing demand of nuclear data of clinically and biologically important 
proton rich radioisotopes of Tc, in this paper, we have investigated two separate routes, e.g., 7 Li 
+ nat Zr and 9 Be + no *Y. Excitation functions of evaporation residues produced in 7 Li+ no *Zr 
and 9 Be -\- nat Y have been measured in the incident energy range 37-45 MeV and 30-48 MeV 
respectively using standard stacked foil technique followed by off line 7-spectrometric studies. 
Other radionuclidic impurities are also produced in these two routes along with Tc. The word 
impurity in the perspective of nuclear data useful in the field of nuclear medicine has been defined 
in our recent paper [32I ]. The impurity describes all the evaporation residues (stable as well as 
radionuclidic) and their decay products including the decay of compound nucleus when aimed 
for a particular radionuclide. Quantification as well as reduction of those impurities by selecting 
proper nuclear reaction parameters like, incident energy, target thickness etc. are important. 
Cross section data provides idea of selecting those parameters. Enriched isotopes can reduce the 
production of impurities, but are quite expensive. Naturally occurring mononuclidic elements are 
preferred as target material to prevent opening of less reaction channels. However, it is also possible 
to use stable elements with more than one naturally occurring isotopes, provided the production 
parameters are controlled. 

Therefore, in the present work, our goal is twofold: 

(i) measurement of excitation functions of evaporation residues produced in 7 Li+™ at Zr and 
9 Be+ na 'Y reactions 

(ii) investigation of the reaction mechanisms involved in those two reactions comparing with the 
nuclear reaction model predictions. 

Broadly, nuclear reaction is understood in terms of three kinds of reaction mechanisms, namely, 
direct (DIR), preequilibrium (PEQ) and equilibrium or evaporation (EQ). DIR reaction is the single 
step process that occurs due to large momentum transfer from projectile to the target nucleons. 
This results in high energy ejectile emission leaving residual at lower excited states. Compound 
nucleus is formed when incoming projectile energy is shared between all the nucleons in the target 
nucleus and equilibrium emissions (EQ) occur due to statistical fluctuation in energy. Between 
these two extremities, PEQ reaction process plays a significant role having features of both. Cross 
section of the product nuclei corresponding to a particular reaction is the resultant of all the 
reaction mechanisms involved. In the present work, contribution from DIR reaction processes is 
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not expected because of low incident energies. We have tried to explain the experimental cross 
sections of evaporation residues in terms of PEQ and EQ reactions in the incident ener gy r ange 
37-45 MeV and 30-48 MeV respectively, using nuclear reaction model codes PACE-II [33j and 
ALICE91 



34, 



35|. 



Section II presents the discussions about the protocol and formalism of the code PACE-II and 
ALICE91 required for the theoretical calculation of excitation functions. Experimental procedure 
is described in section III and section IV deals with the results and discussion of the present work. 



II. NUCLEAR REACTION MODEL CALCULATION 
A. PACE-II 

The code PACE-II is the modified version of the Monte-Carlo code Projection Angular- 
momentum Coupling Evaporation. The deexcitation process of the excited nuclei is calculated 
using the modified version of the code JULIAN, which follows the correctprocedure of angular 



momentum coupling at each stage of deexcitation using Hauser-Feshbach 36| model. The trans- 



mission coefficients for light particle emission of neutrons, protons and as are determined from the 



optical model potential where all the optical model parameters are taken from ref 37|. The shift 
in coulomb barrier during deexcitation is accounted by calculating the transmission coefficients at 
an effective energy determined by the shift. The code internally decides level densities and masses 
it needs during deexcitation. Because of low excitation energy, Gilbert-Cameron level density pre- 
scription is used in the present work with a, level density parameter equals to A/12 MeV^ 1 . The 
ratio of a,f/a n is chosen as unity. Fission is considered as a decay mode. The finite range fission 
barrier of Sierk [38| has been used. Compound nuclear fusion cross section is determined by using 
the Bass method [39j. Yrast parameter is taken as unity. A non-statistical yrast cascade 7-decay 
chain has been artificially incorporated to simulate 7 multiplicity and energy results. 



ALICE91 



The code ALICE91 



34 l35j has been used to calculate the excitation function of product ra- 



dionuclides. Geometry dependent hybrid model 
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42 1 has been used to calculate PEQ 



emissions and Weisskopf-Ewing formalism 43( for EQ emissions. The hybrid model is the combi 



nation of exciton model 44| and Boltzmann master equation approach 45j, |46| ] . It assumes that 
the target-projectile composite system proceeds through two body interaction process. Each stage 
of the relaxation process is designated by the total number (n) of excited particles, i.e., sum of the 
excited particles (p) and holes (h). In each two body interaction, p-h pair may be created or an- 
nihilated or redistribution of energy takes place without changing the number. Hybrid model uses 
never come back approximation, i.e, the model assumes only p-h pair is created in each interaction. 
Hybrid model explicitly determines the pre-emission energy distribution of the excited particles 
which helps to estimate high energy emissions more accurately. Geometry dependent hybrid model 
includes the nuclear surface effects 35J, |4l|, |42J . The PEQ emission cross section for a particular 
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ejectile x with energy e x is given by 



a PEQ (e x ) = £ £ (21 + l) Tl E a D n [/*^SSrJ x^ffej (1) 

Here, A is the de-Broglie wave length of the projectile, 7) is the transmission coefficient of the 
I th partial wave, D n is the depletion factor of the n th exciton state, that is, probability of reaching 
n exciton state without prior emission and f% is number of x type excited nucleon present in 
it. The number hq and n are the initial and equilibrium exciton numbers respectively. The ratio 
N n (l, U, e x ) /N n (l, E c ) is the probability of finding x type nucleon in the n exciton state with energy 
(e x +B x ) where B x is the separation energy of x. The factor A ( e A )^!^( e ) is the emission probability 
of x with energy e x . At(e x ) is the two-body interaction rate. The emission rate X c (e x ) is calculated 
by H 



(2S X + ^H'xtxOinv i^x) 

n 2 h 3 g 

where, S x is the intrinsic spin of x, fi x is the reduced mass, u, M is the inverse cross section of the 
ejectile x with energy e x being absorbed by the residual and g is the single particle level density of 
the composite nucleus. EQ emission cross section is calculated using Weisskopf-Ewing formalism 
as 

&comp T T (3) 

Ccomp is the compound nuclear formation cross section, a is the level density parameter and U 
is the available excitation energy of the compound nucleus after the PEQ emissions. a comp is 
calculated as a comp — a a b s — <Jpeq, where a a b s is the absorption cross section of the projectile in 
the target and ctpeq is the total PEQ emission cross section. 



The calculations have been performed using the code ALICE91 [34|,|35| with geometry dependent 
hybrid model for PEQ emissions and Weisskopf-Ewing formalism for EQ emissions, n, p, d and 
a emissions are considered from the residual nuclides of 12 mass unit wide and 10 charge unit 
deep including the composite nucleus. Fermi gas level density has been used for the calculation 
of reaction cross sections. Reverse channel reaction cross sections have been calculated using the 
optical model. The level density parameter, a is taken as A/12. Rotating finite range fission 
barriers of Sierk has been chosen. Total number of nucleons in the projectile has been chosen as 
the initial exciton number for the PEQ cross section calculation. 

III. EXPERIMENTAL PROCEDURE 
A. Irradiation parameters 

1. 7 Li+ nat Zr 

Pure metallic Zr foil (99.94%) was procured from the Johnson Matthey & Co. Limited. Natural 
assey of Zr contains 90 Zr (51.45%), 91 Zr (11.22%), 92 Zr (17.15%), 94 Zr (17.38%) and 96 Zr (2. 
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To measure the excitation function, self-supporting metallic Zr foils of about 3 mg/cm 2 thickness 
were prepared. Experiments were carried out to measure excitation functions of the product 
radionuclides at BARC-TIFR Pelletron at Mumbai, India. The stack of target-catcher assembly 
was prepared by placing a Zr foil followed by an aluminum catcher foils of thickness 1.5 mg/cm 2 
and was bombarded by 7 Li 3+ and total charge of 1425 /iC was collected over 6.25 h duration. 
Excitation functions of evaporation residues of 93 . 94 > 95 . 96 Tc, 93m Mo and 90 > 96 Nb were measured in 
the projectile energy range 37-45 MeV. 



2. 9 Be+ nat Y 

The only naturally abundant isotope of yttrium is 89 Y, which was procured from Alfa Aesar 
(99.9% pure). Thin self supporting yttrium foils of thickness 3-3.8 mg/cm 2 were prepared to 
measure the excitation function. Target stack was assembled by placing aluminum catcher foils of 
thickness 1.5 mg/cm 2 in between Y foils. Target assembly was bombarded by 9 Be 4+ beam and 
total charge of 304 /zC was collected over 3.75 h duration. Excitation functions of 93 > 94 > 95 Tc and 
93m Mo were measured in the 30-48 MeV incident energy range. 



B. Measurement of activity 

The aim of this paper was to study the production f 93,94 ' 95 ' 96 Tc radionuclides and the associated 
radionuclidic impurities produced from both the reactions. The duration of the irradiation time 
was chosen according to the beam intensity and the half life of the product nuclides. Target and 
aluminum foils were mounted on aluminum ring of 10 mm inner and 22 mm outer diameter with 
0.5 mm thickness. The residual products, if any, recoiled in the beam direction, were completely 
stopped in the aluminum backing. Large area of the catcher foils ensures the complete collection 
of recoiled evaporation residues. At the end of the bombardment, foils were counted for the 7- 
ray activity of the evaporation residues by an HPGe detector of 2.13 keV resolution at 1332 keV 
coupled with a PC based MCA, PCA2 (OXFORD). Efficiency calibration of the detector was 
performed as a function of 7-ray energy using a standard 152 Eu (T 1 / 2 = 13.506 y) source of known 
activity. Each foil was counted for 300 seconds in the live time mode leaving proper cooling time 
after the bombardment and successive measurements were carried out for sufficiently long time in 
the same geometry. In general, projectile energy at a target is the average of incident and outgoing 
beam energy. Beam energy degradation in the target and the catcher foils was calculated using the 



Stopping and Range of Ions in Matter (SRIM) 47J ■ Energy loss is less than 2% in case of zirconium 



targets whereas it is about 4-5% for yttrium targets. Total charge of each irradiation experiment 
was measured by an electron suppressed Faraday cup stationed at the rear of the target assembly 
which in turn is the measure of beam intensity. Evaporation residues recoiled to the catcher foils 
were also measured by 7-ray counting. It was noticed that minute amount of radionuclides recoiled 
to the catcher foil only at higher incident energies and there by neglected in the present calculation. 

The nuclear spectroscopic data of the radionuclides studied in the present work is enlisted in the 
Table U [4^ . The product yields of the evaporation residues in each foil were calculated using the 



G 

I — I. 

reported nuclear data [48]; half life, 7-energy, branching ratio, etc. Background subtracted peak 
area count correspond to a particular 7-ray energy of the 7-spectra is the measure of the yield. The 
product yield (yj) of particular radionuclide (evaporation residue) i at the end of bombardment 
was calculated from the standard relation 



y l = ^ T (4) 

C 7 7 

where c(t) is the count rate at any time t, and are the detection efficiency and branching 
intensity of the characteristic 7-ray of the evaporation residue, designated by i, with decay constant 
X 1 , t is the cooling time. Cross section of the i th evaporation residue (a l (E)) at an incident energy, 
E is calculated from the activation equation 

y* =I p a i (E)n tg x tg (l-e- XiT ) (5) 

where I p is the intensity of the projectile, n tg and x tg are the number of target nuclei per unit 
volume and target thickness, respectively, T is the duration of irradiation. 



C. Uncertainties in measurements 



The expected errors associated in the cross section measurement are as follows: 

(i) Maximum error in efficiency calibration of HPGe detector « 2% 

(ii) Maximum error in determining target thickness (nt g xt g ) in atoms/cm 2 rj 5% 

(iii) Systematic error in the beam current that propagated to the cross section data « 10-12% 

(iv) Uncertainty in the incident beam energy at the successive targets may occur due to the energy 



degradation in the aluminum catchers. According to ref 49, 50], the energy straggling is expected 
to be small even in case of lowest incident energy and hence was neglected in the present work. 

Apart from the above, error occurs in the cross section data due to counting statistics. The 
total associated error related to the cross section measurement was determined considering all the 
factors discussed and the data presented up to 95% confidence level. 



IV. RESULTS AND DISCUSSION 



The residual radionuclides produced through various reaction channels of 7 Li+ na *Zr and 9 Be+Y 
reactions at the maximum incident energies of 44.6 MeV and 47.5 MeV, respectively are shown in 
Figs. [T]and[2]with their characteristic 7-rays. Experimental excitation functions of the residues are 
compared with the theoretical predictions of PACE-II and ALICE91 in Figs. [3j6l Experimental 
cross sections are represented by various symbols with the associated error and theoretical predic- 
tions are shown by the lines (solid line for PACE-II and dashed line for ALICE91). Though the 
code ALICE91 takes care of the PEQ emissions, it has been observed that PEQ reaction has al- 
most no contribution, except about 5% in the highest incident energy, in producing the evaporation 
residues in the energy range studied in the present work. 
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TABLE I: Nuclear spectrometric data of the radionuclides produced through 7 Li+ nat Zr and 9 Be+ 89 Y 
reactions. Bold 7-rays are used in determining the excitation functions. 



Product nuclei 


Spin 


T1/2 


Decay mode (%) 


E 7 (keV) (I 7 %) 


93 Tc 


9/2+ 


2.75 h 


e(100) 


1363.09 (66) 


94m rp^ 


2+ 


52.0 m 


e(100), IT(<0.1) 


871.097(94.2) 


94 Tc 


7+ 


293 m 


e (100) 


7O2.63(99.6),871.097(99.9) 


95 Tc 


9/2+ 


20.0 h 


e(100) 


765.789(94.0), 1073.713(3.75) 




4+ 


51.5 m 


IT (98), e(2) 


778.196(1.9), 1200.165 (1.09) 


96 Tc 


7+ 


4.28 d 


e (100) 


778.196 (99.76), 812.54(82) 


93m Mo 


21/2+ 


6.87 h 


IT(99.88),e(0.12) 


263.143(56.7), 684.753(99.7) 


96 Nb 


6+ 


23.35 h 


P~ (100) 


568.86(56.8) 


90 Nb 


8+ 


14.6 h 


e (100) 


1129.195(92.7) 



are produced through ™ a 'Zr( 7 Li, xn), x being the number of reaction channels. Measured excita- 
tion functions of 93 > 94 Tc and 95 - 96 Tc are compared with theoretical predictions in Figs. [3]and|4] 
respectively in the 37-45 MeV energy range. nat Zr contains five naturally abundant Zr isotopes; 
90 Zr (51.45%), 91 Zr (11.22%) , 92 Zr (17.15%) , 94 Zr (17.38%) and 96 Zr (2.8%). The theoretical 
excitation functions reported in the Figs. [3] - [5] are calculated taking the weighted average of all 
five natural isotopes of Zr. 

It is observed from Fig. [3] that experimental cross sections of 93 Tc are well reproduced by the 
prediction of PACE-II whereas ALICE91 overpredicts throughout the range with a maximum of 
75% at 37.5 MeV incident energy. In case of 94 Tc, both the theoretical calculations are close to 
each other and slightly overpredict the experimental data. However, overall observation shows a 
good agreement between the experimental results and PACE-II predictions. 

Figure [4] shows comparison between the measured cross sections of 95 > 96 Tc and the theoretical 
estimations. In both the cases, experimental data agree with the PACE-II calculations. However, 
ALICE91 overpredicts the data throughout the range. This observation reveals that 93 . 94 > 95 . 96 Tc 
radionuclides are produced as a result of complete fusion of 7 Li in the not Zr target which is according 
to the expectation. It is observed from Figs. [3] and 2] that PACE-II calculations reproduced the 
measured excitation functions with a tendency of slight overprediction, which implies the prominent 
role of compound nuclear reaction mechanism. 

Figure [5] presents the experimental excitation function for the production of 93m Mo and 90 > 96 Nb 
through 7 Li+ na 'Zr reactions along with the theoretical values. PACE-II underpredicts measured 
excitation functions of 93m Mo and 90 Nb by 25% and 50% respectively, at lowest incident energy. 
ALICE91 calculation unexpectedly overpredicts the experimental excitation function of 93m Mo 
by a factor of 2, whereas it agrees with the cross sections of 90 Nb. Underprediction of cross 
section values by PACE-II might be an indication of the PEQ process in formation of 93m Mo and 
90 Nb. However, none of the theoretical models predicts the production of 96 Nb when its signature 
observed experimentally. 
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B. 9 Be+ 89 Y 

Due to the bombardment of 9 Be on 89 Y, mainly 93 > 94 > 95 Tc are produced in the 30-48 MeV incident 
energy range. Only a small amount of 93m Mo is produced at the 47.5 MeV incident energy. Figure 
|6] shows the experimental excitation function of 93 > 94 > 95 Tc along with the theoretical predictions 
of PACE-II and ALICE91. Measured excitation functions of 94 Tc and 95 Tc are well evaluated 
by PACE-II calculations whereas ALICE91 overpredicts the measured cross section of 94 Tc. The 
incident energy range reported (30-48 MeV) in the present work covers maximum production of 
95 Tc and 94 Tc at the lowest and highest incident energy side. The trend of the measured excitation 
function exactly matches with the PACE-II predictions. Measured excitation function of 93 Tc is 
compared only with the ALICE91 values, which starts from 40 MeV incident energy and slightly 
overpredicts the experimental cross sections at next two energy values. However, cross section 
has been measured at 39 MeV, which is well above the threshold value, 36.5 MeV. According to 
PACE-II, 89 Y( 9 Be,5n) reaction channel will open from 45 MeV. It is clear from the cross section 
values that Tc radionuclides are produced through 89 Y( 9 Be, xn), (n=3,4,5), reaction channels 
by complete nuclear fusion. This is expected when the maximum incident energy is less than 
5.5 MeV/nucleon. No signature of PEQ process has been observed from the experimental cross 
sections. Production of 93m Mo is observed only in the first foil at 47.5 MeV. The measured cross 
section is (13 ± 2.1) mb, which was expected from both the theoretical estimations. 

The basic difference in the two theoretical predictions from PACE-II and ALICE91 is observed 
due the formalisms adopted for the simulation of compound nuclear processes. ALICE91 uses faster 
approach of Weisskopf-Ewing sacrificing rigor in physics whereas PACE-II takes Hauser-Feshbach 
formalism sacrificing computing time. The later one found to be better in explaining the nuclear 
data obtained entirely from the compound nuclear reactions. The present work agrees with this 
fact. 

V. CONCLUSION 

Present work reports for the first time two production routes of proton rich Tc radionuclides via 
heavy ion induced reactions of 7 Li+ raa< Zr and 9 Be+ 89 Y. Excitation functions of the evaporation 
residues produced from those reactions have been measured in energy range 37-45 MeV and 30- 
48 MeV respectively. The measured cross section values were found to be in good agreement 
with the detail Hauser-Feshbach model calculations using PACE-II. Cross section data essentially 
revealed the compound nuclear reaction mechanism as per expectation in the energy range studied. 
Experimental nuclear reaction data is important for the quality production of Tc radionuclides. 

In case of 7 Li+ raa *Zr reaction, reported incident energy range, 37-45 MeV, is not sufficient to get a 
complete picture of excitation functions of evaporation residues. More experimental data is needed 
at both, higher as well as lower incident energies. However, presence of all other radionuclides, 
which may cause difficulty in optimizing production parameters of a particular isotope of interest 
reducing other radionuclidic impurity, is prominent from the study. The overall knowledge of 
the work will be helpful in producing clinically important proton rich Tc radionuclides. Cross 
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FIG. 3: Comparison between measured production cross sections of Tc and Tc with theoretical pre- 
dictions of PACE-II and ALICE91. Cross section values are divided by 10 for 94 Tc. 
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FIG. 4: Comparison between measured production cross sections of Tc and Tc with theoretical pre- 
dictions of PACE-II and ALICE91. Cross section values are divided by 10 for 95 Tc. 
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FIG. 5: Comparison between measured production cross sections of 93m Mo and 9C ! 
predictions. Cross section values are multiplied by 0.5 for 96 Nb. 
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FIG. 6: Comparison between measured production cross sections of 95 Tc, 94 Tc and 93 Tc with theoretical 
predictions of PACE-II and ALICE91. Cross section values are divided by 10 for 94 Tc. 
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section data will also help to increase radionuclide purity optimizing nuclear reaction parameters. 
However, it is once again clear that mononuclidic target has advantage in this context. 
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